%% 再入段经纬高转化为XYZ
clc
clear
load('F1.mat')
lla = [lon*pi/180, lat*pi/180, alt*1000];
[SZ, ~] = size(lla);
XYZ1 = nan(3,SZ);
for i = 1:SZ
   XYZ1(:,i) = lla2ECEF(lla(i,:));
end
load('F2.mat')
lla = [lon*pi/180, lat*pi/180, alt*1000];
[SZ, ~] = size(lla);
XYZ2 = nan(3,SZ);
for i = 1:SZ
   XYZ2(:,i) = lla2ECEF(lla(i,:));
end
load('F3.mat')
lla = [lon*pi/180, lat*pi/180, alt*1000];
[SZ, ~] = size(lla);
XYZ3 = nan(3,SZ);
for i = 1:SZ
   XYZ3(:,i) = lla2ECEF(lla(i,:));
end